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We present a detailed study of the statistical properties of a system of diffusing aggregating 
particles in the presence of a steady source of monomers. We emphasise the case of low spatial 
dimensions where strong diffusive fluctuations invalidate the mean-field description provided by 
standard Smoluchowski kinetic theory. The presence of a source of monomers allows the system to 
reach a statistically stationary state at large times. This state is characterised by a constant flux of 
mass directed from small to large masses. It therefore admits a phenomenological description based 
on the assumption of self-similarity and constant mass flux analogous to the Kolmogorov's 1941 
theory of turbulence. Unlike turbulence, the aggregation problem is analytically tractable using 
powerful methods of statistical field theory. We explain in detail how these methods should be 
adapted to study the far-from-equilibrium, flux-dominated states characteristic of turbulent systems. 
We consider multipoint correlation functions of the mass density. By an exact evaluation of the 
scaling exponents for the one and two-point correlation functions, we show that the assumption 
of self-similiarity breaks down at large masses for spatial dimensions, d < 2. We calculate non- 
rigourously the exponents of the higher order correlation functions as an e-expansion where e — 2 — d. 
We show that the mass distribution exhibits non-trivial multiscaling. An analogy can be drawn 
with the case of hydrodynamic turbulence. The physical origin of this multiscaling is traced to 
the presence of strong correlations between particles participating in large mass aggregation events. 
These correlations stem from the recurrence of diffusion processes in d < 2. The analytic methods 
developed here will have more general applicability beyond the study of this specific problem. 

PACS numbers: 47.27.Gs,05.10.Cc,05.70.Ln,68.43. Jk 



I. INTRODUCTION 

The development of a coherent formalism that enables us to calculate the statistical properties of complex systems 
far from equilibrium remains a vexing problem in theoretical physics. It may be that the set of non-equilibrium 
complex systems is too diverse to admit a unified description analogous to the Gibbs formulation for equilibrium 
systems. Nevertheless certain subsets of non-equilibrium systems share enough common features that it is reasonable 
to hope that a general understanding can be obtained. One such subset which has attracted much interest, both 
theoretical and practical, over the years is the class of problems of turbulent type. In the present context, turbulence 
does not refer exclusively to hydrodynamic problems although the statistics of high Reynolds number fluid flow is a 
very important and challenging example. Rather, we take the word turbulent to refer to a class of non-equilibrium 
problems coming from a diverse range of areas in theoretical physics (from hydrodynamics to condensed matter physics 
to aggregation and even cosmology) sharing certain common features which we now describe. 

The defining characteristic of a turbulent system is the existence of a stationary state with widely separated sources 
and sinks of some conserved quantity. While stationary states of turbulent systems lack detailed balance, they are 
characterised by a flux, mediated by nonlinear interactions, of a conserved quantity from source to sink. Universal 
statistics are expected in regions far from sources and sinks, also known as the inertial range. The best known 
example, as mentioned already, is the solution of the Navier-Stokes (N-S) equations at high Reynolds number with 
large scale forcing 1]. It is characterised by an energy flux from large length scales to small length scales where 
viscosity dissipates the energy. Other examples include Burgers turbulence [2| , the Kraichnan model of passive scalar 
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advection 2], kinematic magneto hydrodynamics |2|, and wave turbulence |3(. A common feature of these systems 
is that they each admit a phenomenological description based on constancy of the flux and the assumption of self- 
similarity. The first such theory was the Kolmogorov 1941 (K41) theory of N-S turbulence. Obtaining a quantitative 
understanding of the limitations and applicability of such phenomenology is one of the major theoretical challenges. 
Analytic progress in hydrodynamics has been slow despite strong numerical and experimental evidence for many 
interesting and nontrivial violations of K41-style self-similarity. This is because the N-S equations lack an obvious 
small parameter that would permit a perturbative treatment. Studies of other systems, notably Burgers equation and 
the Kraichnan model, have been more successful 2], providing insight into the breakdown of self-similarity. 

In a recent paper 4] we added to this small class of analytically tractable turbulent systems by showing breakdown 
of self-similarity in a system of diffusing, aggregating masses with a steady source (see section^for definition). In this 
paper we provide a detailed presentation of the results of 4] and the methods developed to analyse this problem. We 
also include several extensions of the ideas. These include the statistics of one-point moments of the mass distribution 
which can be shown to exhibit strong, but in a sense trivial, anomalous scaling and a section which establishes a 
heuristic connection between statistical field theory and an effective kinetic theory description of aggregating particles 
(Smoluchowski approximation). We make extensive use of statistical field theoretic techniques which, while well 
known in the context of equilibrium condensed matter physics, we felt would benefit from some further explanation in 
the context of turbulent systems. The main thrust of the calculations are explained in the main body of the article. 
Several appendices are provided which discuss technical details and aspects specific to the aggregation model which 
we study. 

The layout of the article is as follows. In section |H] we define the "mass model" which underpins our study of 
aggregation. Next, in section lllll we begin where every study of a turbulent system should begin, namely with a 
simple investigation of the implications of dimensional analysis. In this section we explain the K41 phenomenology in 
the context of aggregation and derive the scaling properties for the mass statistics which follow from the assumption 
of self similarity. We then move on to the technical aspects of the article. In section Hvl we give a detailed explanation 
of how the mass model can be described in terms of an effective field theory using Doi's formalism. This can, in 
turn, be mapped into a variant of the A + A — > A reaction diffusion model, a model which has been extensively 
studied. In section we present results from numerical simulations and exact computations of the scaling exponents 
of the first and second order correlation functions. These confirm the breakdown of self-similarity in the mass model. 
The following two sections contain a detailed exploration of this breakdown using renormalisation group arguments. 
Section IVII outlines the structure of the perturbation theory and Feynman rules for the effective field theory and 
outlines how to perform the computations of the relevant diagrams as expansions in e = 2 — d. In section IVIII we 
use reasonably standard RG arguments to derive the Callan-Symanzik equation for the n-point correlation functions. 
We solve these equations in d < 2 to calculate the anomalous scaling curve for the correlation functions as a first 
order epsilon-expansion. We also solve the RG equations in d — 2 to obtain the anomalous logarithmic corrections 
to the mean field answers which one expects to find at the critical dimension. In section lYlIII we return to a more 
physically-motivated discussion and outline the connection between our renormalised perturbation theory and an 
effective kinetic theory description based on the Smoluchowski approximation. We conclude with a brief summary of 
the main results and their importance. Three appendices are provided which supplement the main text with further 
technical explanations necessary for a complete understanding of the work. With the inclusion of the appendices we 
hope that the article becomes completely self contained. 

II. DEFINITION OF THE MODEL 

Consider a cubic lattice in d dimensions whose lattice sites are occupied by particles of integer masses. Multiple 
occupancy of a given site is permitted. Let the number of particles of mass m on site x at time t be denoted as 
A^(x, m). At a given moment of time a configuration of the system is determined by specifying the set of occupation 

numbers, {A^x^, m)}™|^ d . These configurations change in time due to three processes: diffusion, aggregation and 
input. These processes are described below. 

• Diffusion 

JV t (x,ra) -> iVt(x,ra)-l, 
A^(x + n,ra) iV t (x + n,ra) + 1, 

where x + n denotes a nearest neighbour of site x. This rule describes the diffusive hopping of a particle at site 
x to the neighbouring site x + n. The rate at which such hopping occurs is DN t (~x,m)/(2d). Here D is the 
diffusion constant, assumed independent of the mass, and 2d is the coordination number of a cubic lattice in d 
dimensions. 
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• Aggregation 

7V t (x,rai) -> 7V £ (x,mi)-l, 

7V £ (x,m 2 ) -> iV t (x,m 2 )-l, 

]V t (x,mi+m 2 ) — > iV t (x,mi +m 2 ) + 1. 

This rule describes the aggregation of two particles of masses mi and m 2 at the same site x to form a single new 
particle of mass mi + m 2 at site x. The rate at which such aggregation events occur is AiVt(x, mi)iVt(x, m 2 ) 
where A is the aggregation rate, also assumed independent of mass. 

• Injection 

iVt(x, m ) -> iV t (x,m ) + l. 

This rule describes the injection of new particles of mass mo at random sites throughout the system. We take 
mo = 1. The rate of injection of new particles is J. 

We will call the above model the mass model (MM). The MM has 3 physical parameters: the diffusion constant 
the aggregation rate A, and input mass flux J. In addition there are two parameters associated with the lattice - the 
lattice spacing, dx, and the smallest mass, mo- 

Various aspects of this model has been studied in different contexts. The model with parallel update rules was 
initially studied as a simple model for river networks [Sj. A comparison of different physical quantities in this model 
relevant for river networks can be found in 6]. The solution to the single site mass distribution in one dimension 
can be found in (jj, |S( • From the exact solution of the two point correlation function, the exponents governing the 
single site mass distribution was found for arbitrary dimensions The single site mass distribution was solved for 
mass dependent aggregation kernel using Zakharov transformation jlflf and using the Smoluchowski approximation 
in the context of submonolayer epitaxial thin film growth Some other contexts in which the model has been 

studied include granular bead packs , nonequilibrium phase transitions [IH . and self organized criticality . 
All these studies focused on the average mass distribution and its behaviour for large masses. In this paper, we will 
be focussing mainly on multi-point correlation functions. 



III. DIMENSIONAL ANALYSIS AND SELF-SIMILARITY CONJECTURES 
A. Dimensional Considerations for the Mass Spectrum 

Before doing detailed calculations, we begin by asking what we can learn about possible stationary states of the 
mass model from simple dimensional considerations. The first quantity of interest in characterising the long time 
behaviour of the model is the stationary mass spectrum, (N m ). (N m ) is the number of particles of mass m per unit 
volume in the stationary state. Specifically, we would like to know how (N m ) scales with m for large values of m. Since 
we are hoping for universal behaviour in the limit of large m, we assume that the mass spectrum does not depend 
on the position of the source, mo- This assumption must be verified at a later stage. The remaining dimensional 
parameters upon which N m could, in principle, depend are J, D, A and , of course, m. We shall perform most of our 
computations in units where D = 1 but for the purposes of dimensional analysis we shall keep D. 

The dimension of N m is [N m ] = M~ 1 L~ d . The dimensions of the other parameters are: [J] = ML _d T _1 , [D] = 
L 2 T _1 and [A] = L d T _1 . It is immediately evident that there are too many dimensional parameters in the model to 
uniquely determine the mass spectrum. One can readily verify that for any scaling exponent, x, and dimensionless 
constant, Ci, the formula 

(3 — 2x)d (d + 2)x-2d-2 

(N m )=c 1 J x - 1 D^^\ m- x , (1) 

is a dimensionally correct expression for (N m ) which scales as m~ x . This is different to the dimensional argument 
used by Kolmogorov in his theory of hydrodynamic turbulence. For that system, there is a unique dimensionally 
correct combination of parameters giving the energy spectrum. Eq. allows us to pick out the scaling exponent, x, 
for the reaction and diffusion limited regimes: 

* KZ = §, (2) 

K41 _ 2d + 2 ( o\ 
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The above two exponents correspond to a different balance of physical processes in order to realise a stationary state. 
We briefly discuss each to explain the choice of nomenclature. 

• We shall call x KZ the Kolmogorov-Zakharov (K-Z) exponent since it is the analogue for aggregating particles 
of the K-Z spectrum of wave turbulence 0] in the sense that it is obtained as the stationary solution of a mean 
field kinetic equation. This spectrum describes a reaction limited regime where diffusion plays no role. 

• We shall call x K41 the Kolmogorov 41 (K41) exponent since it is a closer analogue of the 5/3 spectrum of 
hydrodynamic turbulence originally proposed by Kolmogorov in his 1941 papers on self-similarity in turbulence. 
This is because in the Navier-Stokes equations there is no dimensional parameter like the reaction rate controlling 
the strength of the nonlinear interactions. This exponent describes a diffusion limited regime where the reaction 
rate, A, plays no role, reactions being effectively instantaneous. 

The case x = 1 corresponds to one in which (N m ) does not depend on the mass flux J. However, this is not of 
physical interest for this problem. On the other hand, each of the regimes characterised by x KZ and x K41 carry a 
mass flux and is relevant. 



B. Self-Similarity Conjectures and Multipoint Correlation Functions 

We are interested in more than just the average mass density (N m ). To characterise correlations in the the mass 
model we must also consider multipoint structure functions. Let C n (rai, . . . , m n )(AV) n Y\ i drrii be the probability of 
having particles of masses rrii in the intervals [m^m^ + drrii] in a volume AV for i = 1 . . . n. C\(m) is the average 
mass density (N m ). We ask how C n (mi, . . . , m n ) varies with mass when mi, . . . , m n ^> mo- In particular what is the 
value of the homogeneity exponent j n defined through C n (Tm\, . . . , Ym n ) = T~ ln C n (mi 1 . . . , ra n )? 

As for the density, dimensional analysis alone is insufficient. The formula 

(3n-2 7Tt )d (d + 2^ n + (2d + 2)n 

C n (mi,...,m n ) = c n J ln n D d-2 \ d-2 

(mi . . .m n )~°£, (4) 

is dimensionally consistent for any exponent 7 n where c n is a dimensionless constant. We are again assuming that 
the large mass behaviour of the C n 's is independent of tuq. The simplest way to obtain a theoretical prediction for 
the mass dependence of the C n 's is to use a self-similarity conjecture similar to Kolmogorov's 1941 conjecture about 
the statistics of velocity increments in hydrodynamic turbulence. Assume that C n depends only on the masses ra$, 
mass flux J and either the reaction rate, A or the diffusion coefficient D. Depending on which assumption we make, 
dimensional analysis allows us to determine the mass dimension of C n . For the reaction limited case we obtain 

7„ KZ = §», (5) 



and for the diffusion limited case, 

K41 2d + 2 

/n d + 2 w 

Note that in both cases, the dependence of 7 n on the index n is linear, reflecting the assumed self-similarity of the 
statistics of the local mass distribution. When n = 1, 7f" 41 = (2d + 2)/(d + 2) coincides with the result of an exact 
computation of 71 for d < 2 [3 so we expect that the K41 conjecture is the appropriate theory in d < 2. In d > 2 
it is known that 71 = 3/2, hence 7^ z = 3/2 is the correct scaling for the density. Therefore the KZ conjecture is 
appropriate in higher dimensions. For d > 2, the statistics of the MM should be accurately described by mean field 
theory (KZ) and the self-similarity conjecture should hold. In this paper we will concern ourselves with d < 2. 

The K41 self-similarity conjecture assumes that C n does not depend on A, mo, the lattice spacing, and the box 
size AVdmi . . .dm n . The lack of dependence on the lattice spacing is expected due to the renormalizability of the 
effective field theory describing the MM below two dimensions. We will however find an anomalous dependence of 
correlation functions on a length scale depending on the other parameters and the box size which leads to a violation 
of self-similarity. 



5 



IV. FIELD THEORETIC DESCRIPTION OF THE MASS MODEL 



A. Motivation 



In order to check the validity of Kolmogorov conjecture for the mass model for n > 1, we need to go beyond 
dimensional analysis. We shall do this by constructing an effective field theory that provides us with a continuum 
description of the model. We are helped by the fact that the effective field theory which describes the mass model is 
closely related to a much simpler and well studied theory which describes the A + A — > A reaction-diffusion model. We 
shall then use standard techniques of statistical field theory to extract information about the mass model from this 
theory. As is well known, this model has critical dimension 2. In dimensions greater than 2, a mean field description, 
characterised by K-Z scaling works. In d < 2, this mean field description breaks down and correlations between 
particles become important. We shall show how to take into account these correlations and we demonstrate how they 
lead to the breakdown of self-similarity in d < 2 by calculating the scaling behaviour of the correlation functions, 
C n , for large masses. Unlike for the case hydrodynamic turbulence, such an analysis is possible for the mass model. 
This is because in d < 2, the large mass statistics of the model in are governed by a perturbative fixed point of the 
renormalization group flow in the space of coupling constants of the model. The order of this fixed point is e = 2 — 
which allows one to compute the relevant scaling exponents in the form of an e-expansion. 



B. Effective Action for the Mass Model 

Using Doi's formalism, it is possible to construct an effective field theory of the mass model. The steps in the 
procedure are as follows: 

1. Write a master equation for time evolution of P({A^(x^, m)}), the probability of finding the system in a given 
configuration, {A^(x^,m)}. The master equation is linear and first order in time. 

2. Introduce creation and annihilation operators, a^ m and a\ m , which create and destroy particles of mass m at 
site x^. Then convert the master equation into a Schroedinger equation : 

j t m)) = -HWi,m,a{ m ] 

using Doi's formalism (second quantisation) [Tfil lr^. 

3. Use the Feynman trick to derive a functional integral measure which converts the second quantised Schroedinger 
equation into a continuous field theory. 

These steps are well described in the context of reaction-diffusion models in For the model of interest here, 

the procedure is similar to the reaction-diffusion case with the algebra made slightly more complicated by the necessity 
to keep track of a mass index for each particle. Explicit formulae for the master equation and Hamiltonian operator 
of the MM are given in appendix 1X1 

After going over to a path integral formulation of the Schroedinger equation we can express the average density 
and other correlation functions as path integrals. Further detail on this procedure can be found in appendix |S1 For 
example the average density, N t (m) = (</> x ,m(T)), is given, in the notation of appendix iBl by 



N t (m) = J W^mWe-^^^^ 1 , (7) 

where 

Smm[0,0*,A^,A] = J^dr J d d xdm {</>*d t <l> 



+#[</>,</>*]}, (8) 



and 



H[4>,4>*} = DV x ct>* m -V x ct> m + -6(m-m )tt>* m 

m 



A J drriidrri2 {£( m 2 — m — mi) 

[Ci 2 - 2 </4 - OmJ 0m0mi} • 
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C. Dimensional Analysis of the Effective Action 

As for the case of stochastic aggregation without source [2(| , it is helpful to nondimensionalise the fields <p and (jf 
and express the action in eq. (|SJ) solely in terms of dimensionless quantities. Introduce dimensionless fields, </> and </>*, 
in eq. {7J) by the following rescalings: 

r —> t r, 
x — > a/I^x, 
m — > XJt 2 m, 
1 — 



X 2 Jt 3 r ' 

</>* - 0*, 

to obtain 

N t (m) = ypW^^e-^™^*' 1 ' 1 ' 1 ' 11 , (9) 
where the dimensionless interaction coefficient is 

g = Dit^X. (10) 

The fact that g — > as £ — > oo for dimensions greater than 2 and g — > oo as t — > oo for dimensions less than 2 expresses 
the well known fact that the critical dimension of the mass model is 2. We shall have much more to say about this 
later. 

D. The Stochastic Smoluchowski Equation 

It is possible to establish an exact map between the field theory in eq. J7J) and the following stochastic integro- 
differential equation, 0, [2(J : 

g \ r rn 

— DV 2 ) 4>{m) = A / dm! 4>(m!)(j)(m — rri ) 

-2\<j>(m)N + —5(m - m ) + iy/2\(j>(m)r]{x , t), (11) 
m 

where TV = Jq 00 dm f (j)(m f ), z = V^—l? an d w(%)t) is white noise in space and time: 

(rj(x , f)^, 0) = - 0^ - ( 12 ) 

The technical details of this mapping can be found in appendix El Without the noise term, one recognises eq. (jTTJl 
as the mean field (Smoluchowski) equation of the model. Thus, all fluctuation effects are encoded in the imaginary 
multiplicative noise term. 

E. Correspondence with A + A — > A Model 

Eq. simplifies after taking Laplace transform with respect to the mass variable poj . Let 



Then, 



Rn(x,t) = / dm<j)(x,rn,t) — I dm<fi(x,m,t)e /im . (13) 

(J^-W 2 ^(£,£) = -XRl + ^ + 2iy/XR li (S,t)ri(S,t), (14) 

where = J(l — e _/im °). In terms of field t), eq. ((14)) becomes a stochastic version of the rate equation for the 

A + A — > A reaction in the presence of a source. Hence, the computation of the average mass distribution in the mass 
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FIG. 1: Propagators and vertices of the theory. 



model reduces to solving a one-species particle system with a /i-dependent source and then computing the inverse 
Laplace transform with respect to \i. For example, to compute the average density, (iV m (t)), for the mass model, we 
first calculate (i?^), the average of the solution of eq. ([Tl)) with respect to the noise, rj(x,t). We then take the inverse 
Laplace transform with respect to \i and obtain the density from eq. ([T3|). By applying the Martin-Siggia-Rose 
(MSR) procedure 21] to eq. we can write (R^) as a functional integral : 

(R^) = J VR^VR^ R^e~ s ^ R ^\ (15) 

where the effective action for the reaction-diffusion system described by eq. (|14|) is 

Sbd[R»,R»] = Jdxdt R fI (d t -DA)R fI + \R fI Rl 



-\-\R 2 R 2 — —Ru 



(16) 



In order to compute higher order correlation functions C n (m, t), we need to know correlation functions of the form 
(R fI1 (x,t)R fl2 (x,t) . . .R^ n (x,t)). These are non-trivial, as the stochastic fields R^(x^ t)'s are correlated for different 
values of \i via the common noise term in eq. El To clarify what is meant by this, we apply the MSR procedure to two 
copies of eq. (2| describing the evolution of R fI1 and R^ 2 respectively to obtain a functional integral representation 
for (R^R^). This gives : 



(R^R^) — J VR^VR^VR^VR^ R^R^ 



The point to note here is that the path integral measure for correlations of for different values of \i does not factorise 
owing to the presence of the last term in eq. (|17|). Thus, to compute n-point correlation functions in the MM, one 
needs to analyse a system of n stochastic rate equations of A + A — > A theory coupled via common noise terms of the 
form shown in eq. ([T7)). Some detailed explanation of the physical interpretation of higher order correlation functions 
in terms of the probability of multi-particle configurations in the particle system is provided in appendix [0 



F. Feynman Rules 

Solving eq. (|14|) perturbatively in A and j, and then averaging over noise, one can derive the set of Feynman rules 
for the computation of correlation functions. Alternatively they can be written down directly from the action, eq. 
([Trjj) . See pj| for the details of the procedure. However care must be taken to include the "extra" vertex which arises 
when computing correlations between fields with different fi indices. The Feynman rules are summarised in fig. ^with 
time increasing from right to left. The slightly more complicated prefactor for the quartic vertex takes into account 
the aforementioned "extra" vertex. 

The propagator is just the regular diffusive Green's function which, in d spatial dimensions, is 

Go(x 2 -xi,*2-*i) =4tt(*2-*i)~* e"^™, (18) 
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FIG. 2: 7 n as a function of n is one dimension. The straight line shows the Kolmogorov answer [eq. jj]. The dotted line 
shows eq. (|74J) with e = 1 and terms of order e 2 and higher set to zero. The values 70, 71 and 72 are exact. 73 and 74 were 
obtained by Monte Carlo simulations performed on a lattice of size 10 5 and averaged over 2 x 10 7 Monte Carlo time steps with 
J = 4D. 



or 

G„(M) = (2tt)-* e~ k2t , (19) 

in the momentum-time representation usually used in computations. 

The one-point function, (R^) is then given by the sum of all diagrams constructed from the building blocks shown 
in fig. [l]with a single outgoing line. Likewise, the n-point correlation function (R fll (xi,t\) R^ 2 (#2, h) ... R^ n (x n , t n )) 
is given by the sum of the contributions of all diagrams which have n-outgoing lines. In section fVll we shall turn to 
actual computations. 



V. NUMERICAL SIMULATIONS AND MULTISCALING 
A. Numerical Simulations of the Mass Model in d=l 



We first look at the results of Monte Carlo simulations of the MM which confirm that there is indeed some interesting 
behaviour which requires explanation. In particular, numerical simulations show a breakdown of self-similarity in the 
mass model in one dimension and multiscaling of the correlation functions, C n . The results are shown in fig. [21 



B. Constant Flux Relation - Analytic Confirmation of Multiscaling 

We know that the K41 hypothesis works for n = 1. From fig. it is clear that the Kolmogorov scaling breaks down 
for n > 1. It is also possible to analytically confirm that j n ^ 7^ 41 in d < 2 by computing 72. From the definition of 
72, if follows that 

$ 2 ( mi ,m 2 ) = f-L_Y 2/2 J^-), (20) 
\rri1m2 J \rri2 J 

where <p is an unknown scaling function which satisfies <j>(x) = (j)(l/x) due to a symmetry. Our aim is to compute 72 
without using the e-expansion which we shall use in section IVTTI to compute j n for general n. 
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FIG. 3: Diagrammatic form of mean field equations for R^(t) and the tree-level propogator, G™ f (2; 1). 

As we are interested in 3>2(^i> 777,2) for mi, m2 > 0, 

$ 2 (m 1 ,m 2 )=/ / dfud^iR^M^-^e-" 1 ^, (21) 

J a — ioo J a — ioo 

where solves eq. ([T3j). Due to eq. (J20J), 

(^ 1 ^ 2 )=(-V) 1 72/ V(^), (22) 

where -0 is an unknown scaling function. To find the large mi, m2 asymptotics of <£, we need to know the small /ii, /i2 
asymptotics of (i? /Ltl i? /X2 ). Averaging eq. (|T^|). with respect to noise and setting d t {R IJL ) = in the large time limit, 
we find that (R^R^) = « ^ for /i <C mo- Comparing this result with eq. ([22]) we find that 72 = 3. 

Note that 72 does not depend on dimension, <i, of the lattice. Therefore, it is correctly predicted by mean field 
theory. The non-renormalization of 72 by diffusive fluctuations can be explained by mass conservation or, more 
precisely by constancy of the average flux of mass in the mass space, see [lOj for more details. Here, we simply wish 
to point out that the exact answers for 71 and 72 establish multiscaling non-perturbatively: the points (0,0), (l,7i) 
and (2,72) do not lie on the same straight line. 

Due to its close connection with mass conservation, the law 72 = 3 is a counterpart of the 4/5 law of Navier-Stokes 
turbulence. Recall, that 4/5 law states that the third order longitudinal structure function of the velocity field scales 
in the inertial range as the first power of the separation. It is interesting to notice, that Kolmogorov theory respects 
4/5 law in Navier-Stokes turbulence, but violates 72 = 3 in the MM. 



VI. PERTURB ATIVE EXPANSION FOR CORRELATION FUNCTIONS 



A. Mean Field Analysis 



The mean field theory associated with the field theory described by the effective action, eq. © can be thought of in 
several complementary ways. Let us suppose that the reaction rate, A, is the smallest parameter in the problem. This 
means that the dimensionless interaction coefficient, g, given by eq. JTJJJ) is small. In this case, the path integral in 
eq. (0 can computed in the limit g — > using the saddle point method. In this limit, <p satisfies the Euler-Lagrange 
equation (expressed in dimensional variables) : 

/ g \ r rn 

^— — DV 2 j (j){m) = A J dm (j){rn')(j){rn — rn ) 

-2\Mm)N + —S(m - m ), (23) 

which we recognise as the mean field equation derived for classical aggregation problems by Smoluchowski. Now, 
if ^ — > then it follows that the noise term disappears from the non-dimensionalised version of the stochastic 
Smoluchowski equation, eq. JTTJ, leaving us with a deterministic equation for the density, which is again the classical 
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Smoluchowski equation. For readers interested in the analogy between stochastic aggregation and wave turbulence, 
the mean field Smoluchowski equation is the analogue of the kinetic equation. 

In [HI we studied in great detail the stationary state of eq. ([23)) and showed that the spectrum 



jVm = Vi^A m " 1 ' (24) 

is the exact stationary solution as t — > oo. This solution carries a constant flux, J, of mass from small masses to large. 
This is the Kolmogorov- Zakharov spectrum of the mass model which we identified from dimensional considerations 
in section ITTT1 as corresponding to a reaction limited regime. 

When do we expect the mean field answers to be correct? The K-Z solution is established in the limit of large 
times. Since the mean field results become exact in the limit g — > 0, eq. (|10|) implies that the spectrum ([21)1 should 
be correct for d > 2. For d < 2 the mean field approximation quickly breaks down and we must take into account the 
effect of fluctuations. This will be the main objective of the rest of this article. 

Let us now identify clearly the terms in the diagrammatic expansion which give the mean field answers so that 
we can see how to use our formalism to compute the fluctuations about the mean field. Since the mean field kinetic 
equation corresponds to the deterministic limit of the stochastic Smoluchowski equation, the corresponding field 
theory has no loops. Therefore we expect the mean field answers for the average density to correspond to the sum of 
all tree diagrams with a single outgoing line. Let us now analyse these. 

Let Rmf, denoted by a thick line with a cross, be the contribution to R from all tree level diagrams. The equation 
satisfied by R m f is shown is diagrammatic form in fig. In equation form, it reads 



™gL = h-- XRl f . (25) 
at mo 



This is easily solved to give 



R »<<<> =^'™mV™ ( -v^a- (26) 



Performing the inverse Laplace transform in the limit \i — > we find that as m — > oo, 



and recover the K-Z spectrum as we should. Both the constant and the exponent agree with those obtained by the 
Zakharov transformation of the mean field kinetic equation confirming that our approach makes sense. 

It is convenient to define G]^ f (£2^2; xiti) as the propagator that includes all the tree level diagrams. The equation 
obeyed by it is shown in fig. [5f>. The solution is 



Gf(2;l) = G (2;l) 



r— -1 2 

cosh J^h 

V mo 1 



(27) 



tl - 2 ^°° G (2;l)e- Q ^-^, (28) 
where Go is the Green's function of the linear diffusion equation, eq. ([T5)l . and 



fi M = 2 A /^, (29) 
w mo 



is the inverse of the mean field response time. 



B. Loop Expansion 

In order to take into account fluctuations about the mean field answer we need to compute diagrams with loops. 
Ordering the terms in the perturbation series according to the number of loops is known as a loop expansion. Using the 
mean field density and response functions computed above simplifies the task of computing the sum of all diagrams 
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FIG. 4: Zeroth and first order terms in the loop expansions for (R^) and (R^R^) 



with a given number of loops. We now demonstrate by power counting that loop expansion of the mean mass 
distribution corresponds to weak coupling expansion with respect to A. The quantity (R(x,t)) is given by the sum 
of all diagrams with one outgoing line built out of blocks shown in fig. Consider such a diagram containing L 
loops, V vertices and N R m f-\mes. The corresponding Feynman integral contains (in the mixed momentum-time 
representation) dL momentum integrals and V time integrals. Hence, the integration over all times and momenta 
produces the factor n~ v + dL / 2 ~ + . The N i? m /-lines produce the factor R^ f ~ A _7VA/2 . A factor X v comes 

from V vertices of the graph. Hence the corresponding Feynman integral is proportional to + . Note also 
that the number of triangular vertices in the graph is equal to N — 1 and the number of quadratic vertices is equal 
to the number of loops L. Thus the total number of vertices is given by V = L + N — 1. Therefore, any L-loop 
graph contributing to the average mass distribution is proportional to A - 2 + 2 _ ( 1 +2). We conclude that loop expansion 
corresponds to the perturbative expansion of R around the mean field value with the parameter A - * - . 



C. Breakdown of Loop Expansion 

The conditions under which the loop corrections to the mean field answer can be neglected are most simply derived 
using dimensional analysis. The scale of diffusive fluctuations is given by the only constant of dimension length which 
can be constructed out of \i and J: Ld = (/iJ) _1// ^ d+2 ^. The dimensionless expansion parameter in the loop expansion 
above is go(fi) = AL^ 2 , where e = 2 — d. The large mass behaviour of N m is determined by the small- ji behaviour of 
In d < 2, go goes to infinity in the limit /i — > and the loop expansion breaks down. Thus a re-summation of 
the loop expansion is needed in order to extract the large-m behaviour of iV m in low dimensions. 



D. Calculation of One Loop Corrections to Mean Field Theory 



Let us now compute (R^) and {R^R^) to one loop order. The diagrams are shown in fig. 01 The corresponding 
algebraic expressions are evaluated by dimensional regular isat ion in dimension no longer necessarily an integer. 



{R, 



3n 

Amo 



1 + 



(47T) 



(I) 



Amo 



(30) 



where . . . represent terms of higher order in A which necessarily have more loops. In this formula we have introduced 
the quantity e, defined as 



e = 2-d, 



(31) 
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to measure the deviation of the dimension of the system from the critical dimension. 

The four diagrams for (R^R^), shown in fig. 21 give the following respective contributions : 



/ p p \ _ pmf pmf _ 



2A « ^ ^ 



+i? (1) i? mf + i? mf i? (1) 



(32) 



1 



(I) 



an expression which is correct to one loop order. Note that the second diagram in the expression for (R^R^) 
describes the correlation between R^ fields for different values of \i and prevents the factorisation of the 2-point 
function into a product of 1-point functions. We shall need these expressions again to when we use RG to resum the 
loop expansion. 

For a given mass scale, m, there is a corresponding \i scale, 1/m and a corresponding length scale, L M defined as 



( 3jm 
\Dm 



d + 2 



At this point, let us also define the dimensionless reaction rate, g, as 

9 = AL^. 

In what follows, it shall be convenient to express eq. (|3()|) and eq. (|32|) in terms of L M and g. This gives 



L- 2 1 



4(2tt) 1 



(1) 



and 



(-R/Xl^/i 2 ) 




(33) 



(34) 



(35) 



(36) 



Note that the fi dependence of this expression is illusory since g also depends on ji. To study the behaviour of the 
corrections to mean field answers which we have just calculated, we need to study the large m behaviour of the Laplace 
transforms with respect to the /i's of the expressions in eq. (|5U|) and eq. Simple calculation shows that the 

second terms inside the square brackets in these expressions diverge as the /i's are taken to zero when e > signifying 
a breakdown of the loop expansion. This is as expected from the power counting argument of section IVI CI 



VII. RENORMALISATION GROUP ANALYSIS FOR d < 2 



A. Epsilon Expansion 



The loop expansions for correlation functions computed in section lVI Dl are expansions in powers of the dimensionless 
reaction rate, g. The problem is that in d < 2, these expansions become badly ordered as we approach L M — > oo. 
However, since we have computed correlation functions in arbitrary dimension, we can convert the loop expansions 
into expansions in e = 2 — d at a fixed value of g. For e <C 1, eq. (|35|) and eq. can be written : 



L- 2 — 



9 

47T6 



and 



9 

27T6 



(37) 



(38) 
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FIG. 5: Diagrams contributing to the renormalisation of the reaction rate. 



Of course these series are still badly ordered as e — > 0. The idea is to replace certain correlation functions with 
appropriate renormalised quantities, also expressed as expansions in e, such that the renormalised counterparts of 
the above expressions are well-ordered in e. The final pay-off comes when we find that these expressions remain well 
ordered in e even when we take the limit L M — > oo because of the presence of a perturbative fixed point, a structural 
feature of the theory which we must now explain in order to make sense of this scheme. 



B. Renormalised Reaction Rate and ^-function 

The presence of a perturbative fixed point for the A + A — > A model was originally pointed out by Peliti [22^ . 
The corresponding calculations in the presence of a source were done by Droz 23]. This is sufficient to deal with the 
problem at hand. Nevertheless we shall paraphrase their arguments here for the sake of completeness. 

Let us define a renormalised reaction rate, A#, as the amputated 3-point vertex function shown diagrammatically 
in fig. After performing the algebra we find 



2(2tt) 



* + 



(39) 



Now we introduce a dimensionless renormalised reaction rate, gR = XrL 6 , as we did in eq. which is given by 



9R = 9 ~ 

For small values of e this can be written as 

9r 



1 T (i\ 

2(2 7 r) 1 -f V4/ 



where 



2(2tt) 
1 

27T6 



(40) 



(41) 



(42) 



o(l) e < 1. 



Inverting eq. pi*]) allows us to convert perturbative expansions in the bare reaction rate, g, into expansions in the 
renormalised reaction rate, gn- We find 



9 = 9R 



-\e)gj 



(43) 



The crucial point to all of this analysis is the following observation. Although for positive e, g diverges as L M — > oo, 
rendering perturbative expansions in g useless for capturing the large mass behaviour of the theory, we will find that 
gn remains finite as \i — > oo. Furthermore, gR tends to a value which is of order e. Therefore for small e we can use 
eq. (Jl3|) to convert expansions in g into expansions in gR which then have a better chance of remaining non-singular 
when we take L M — > oo. The process of replacing g with gR is usually called coupling constant renormalisation in the 
literature. 

The large mass behaviour of the renormalised reaction rate is determined by the /^-function of the theory defined 



P(9R) 



L 



9gn 
"dL n 



(44) 
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Using the fact that L^gj^g 71 = neg n together with eq. and eq. (|35J) we quickly find 

f3(g R )=eg R (l-g; 1 (e)g R + ...). (45) 

Eq. (HH) now tells us how qr changes as we vary L M . Solving this differential equation with the initial condition 
g(Lo) = go determines how the reaction rate varies with scale. The behaviour is different in d = 2 and d < 2. In 

d < 2, 

so that 

9R( L ») = 7] ffom", So re ' ( 47 ) 

We note that gR goes to a fixed point value of 2tt€ (^corrections of 0(e 2 )) as L M — > oo irrespective of the initial values 
of #o and L . For simplicity we can take L = (1 — giving 

9*W = fffi . (48) 

This universal behaviour as L M — > oo is what is meant when we say that the renormalisation group flow has a 
perturbative fixed point in d < 2. 



In d = 2, e = so that 



which gives 



g R (LJ = g ° , , . (50) 



1 



Thus in d = 2 the reaction rate decays to as L M — > oo but logarithmically slowly and, unlike in the case d < 2, 
retains some memory of the small scale cut-off, Lq. 



C. Average density in d < 2 

Let us now show how all of this technology works by calculating the large mass behaviour of the average density. 
We define the renormalised density, (^)r, by using eq. to replace g with the renormalised reaction rate, gR in 
eq. (|37j). Using eq. (|42|) . a Taylor expansion shows that the replacement of g with gR cancels the e-singular term so 
we are left with 

{R fl ) K = L^^=[l + o{g R )}. (51) 
VSR 



Now as L M — > oo, — > so we can now take the limit to obtain 



<i2 M ) R ~L«- 2 --L[l + ( e 2 )]. (52) 



Finally note that as \i — > 0, 

J/i 



/i 



m 



allowing us to perform the inverse Laplace Transform required to return to mass space. Using the definition, eq. JT3J), 
of we finally find 
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giving a scaling exponent which we know to be correct 9]. We recognise this as the K41 exponent. As discussed 
in section lllll we could have obtained this answer simply from dimensional arguments once we recognised that the 
reaction rate is renormalised away to infinity and hence cannot play any role in the answer. However this would ignore 
the possibility of anomalous dimensions. By calculating the one loop corrections to (R^) we confirmed the absence 
of any relevant (for (R^)) couplings other than the reaction rate itself. We shall find that this is not the case for the 
higher order correlation functions. 



D. Higher order moments of the density 



A natural object to study to gather more information about the mass distribution function would be moments of 
the density of the form M n (m) = (N^). As explained in the appendix 151 ( see eq. IB 181 and the explanation thereafter) 
these moments exhibit "extreme" anomalous scaling characterised by Burgers-like scalings : 

(N») ~ (N m ) mZoo m- 2 -*W. (54) 

For the MM however, this anomaly is somewhat trivial from a physical perspective. It arises because large masses 
become large by absorbing almost all nearby particles. Thus asymptotically, the number of heavy particles on a given 
lattice site ends up being either zero or one. Taking moments of such a distribution will always give the behaviour 
described by eq. ((51)). However, the analysis of appendix |B] which allows one to extract this essentially non-mean 
field behaviour from an initially weakly coupled theory is not trivial and can be expected to yield interesting results 
in other contexts. To observe true the multiscale structure of the mass model one should really study multipoint 
correlation functions. We do this next. 



E. Higher order multi-point correlation functions in d < 2 

The analysis for the higher order correlation functions is not quite so simple as for the density. By replacing g with 
qr in eq. (|36j) we get 



(RuiR^) 



(^i)r(^ 2 )r 



1 



9R , / 2 



R) 



(55) 



We see that we have removed the e-singularities from the (R/j,) factors but the singularity inside the square brackets 
remains. The correct definition of the renormalised 2-point function must include renormalisation of the amplitude 
of C2, not just the reaction rate. This process is known as composite operator renormalisation. The correct definition 
of the renormalised 2-point function is therefore 



{R^R^r — Z2(R^ 1 R l 



(56) 



where the amplitude Z 2 is chosen so that (R^R^r is nonsingular in e: 



Z 2 = 1 



9r 

27T6 



(57) 



The prefactor, Z n , of the n th order correlation function can be computed in a similar manner to the second order 
one. For example, in the loop expansion of the 3 rd order correlation function, there are three diagrams containing 
singularities which are not removed by coupling constant renormalisation. These are shown in fig. El For the n-point 
function there are ^n(n — 1) such diagrams. Each of these diagrams contributes to the one-loop expression for 
Z„ so that : 



Z n = l + -n(n-l) 



9r 

27T6 



o(g 2 R )- 



(58) 



This situation is a bit more complicated than before. To extract the scaling exponent we employ the technology 
of renormalisation group (RG) which was not truly necessary to compute the scaling of the density. Our discussion 
follows closely the presentation of p3 |. The approach is based on the simple observation, already made at the end 
of section IVI Dl that the n th order correlation function, C < ' n \L fI1 . . .L^ n ) = (i? /Xl . . .R^ n ) does not depend on the 
arbitrary length scale L M , known in RG language as the reference scale. It immediately follows that 



d 
dL, 



(z-\g R )C^\L 



0. 
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The L^- dependence of the bracketed expression comes from three sources : an explicit dependence of on L M , an 
implicit dependence through g^L^) and an implicit dependence through Z n (^(L /i )). We can thus write 



- JL. 



dg R d dZ n d 

J-'ii.^z — ^ r L, 



z-\g R ) ) ( V ---Wn, 9R, = 0, 



(59) 



where the partial derivative with respect to is now taken at fixed g R and Z n whose dependences on are catered 
for by the additional derivatives. This can then be arranged to give the equation : 



d 
dL b 



P(9R) 



d 
dg R 



lni.9R) 



(Lm ■ ■ ■ L^ n ,gR,Lfj,) — 0, 



(60) 



where 



L 



»dL L 



(logZ 2 (^fl)) , 



OR 

2tt 



■°(9r), 



(61) 



and P(gR) is given by eq. By itself, this equation just tells us how varies with physically meaningless 

reference scale, L M . However dimensional analysis provides extra information. Since the physical dimension of C ] 
is L~ nd it must satisfy an Euler equation [24j | 



(n) 



nd 



C^ ) (L fll ...L fln ,g R ,L fl ) = 0. 



(62) 



Suppose we now rescale all lengths by some amount, A, by introducing L fIi = AL Mi . Eq. (|62|) allows us to convert 
derivatives with respect to L M into derivatives with respect to A: 



d 



so that eq. (|6()|) can be written as 



~ A dX + ^ 9R ^dgR~ ~ nd ~ ln ( 9R ) 



C^ ) {L IH ...L ltn ,g R ,L li )=0. 



(63) 



(64) 



This equation is called the Callan-Symanzic (C-S) equation. It tells us something physically useful, namely how the 
renormalised correlation function changes as we rescale its arguments by A. We wish to solve it in the limit of large 
A. This can be done using the method of characteristics. For A = 1, is given by the mean-field answer which is 
valid for small values of the L^'s, thus providing an initial condition : 



4 n) (A = 1, g R = g ) = g 3 (L Ml . . . L tln )~ d . 
For d < 2, (3(gR) = egn(l — |^-) and the characteristic equations are 



(65) 



dA _ 
ds 

dgR n gR v 

^7 = egR{1 ~ 2^ } ' 



da 



(n) 

R _ 



ds 



_ 1 , x g 
nd H — n(n — 1) — 
2 v 7 2tt 



(66) 
(67) 

(68) 



with the boundary conditions 



A(# ,s ) = !> 

gR(go,s ) = g , 

C£\g ,s ) = g-^{L^...L^)- d . 



(69) 
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FIG. 6: Diagrams contributing to the anomalous dimension of the 3-point function (R fJ , 1 R^ 2 R fJ , 3 ) at the one loop level. 



If we solve these equations, use the uniqueness of the characteristic curves to express so and go in terms of A and 
gu and then evaluate the solution at s = 0, the answer can be found explicitly: 

n / ( 2 \ x ^ n ( n— 1) 

(^H^lH) 1 ■ ,70> 

Taking A — > oo we conclude 



4»>(L M1 . . . L,„, te> L„) ~ J] 4 B(B_1) V ^ LT d -^ n -'\ (71) 



independent of the value of gR. This independence is the consequence of the presence of a fixed point of the /3-function. 
All values of g^ flow to the fixed-point value, g* = 27re, leaving a universal answer in the limit of large A. It remains to 
perform the inverse Laplace Transform to find the scaling properties of the original mass- space correlation functions. 
To do this we note from eq. that for large values of the Z^, 

It is then easy to perform the n inverse Laplace Transforms with respect to the fli to get 

C ( £\m 1 ...m n ,g R ,L ll ) ~f[m:^"^. (73) 

1=1 

The mass scaling of is therefore m~ lri with 

2d + 2 n(n — l)e , 9s ,^ AS 

In = n— h -Vl T- + )• ( 74 ) 

/n d + 2 2(d + 2) V ; V ; 

Note that j n acquires a correction to the value predicted from K41 theory signalling the breakdown of self-similarity 
in low dimensions. This is the multiscaling curve against which we compared our numerical results in fig. [21 

F. Logarithmic Corrections in d = 2 

In d = 2 scale invariance is broken by the presence of logarithmic corrections to the mean field scaling. For 

completeness, let us calculate the powers of the logarithms acquired by the s. In d = 2, P(gFt) = — §^ and the 
C-S equation, eq. reads : 



-A— - ^- 2d - -n n - 1 — 

<9A 27rG%? 2 v J 2ir 



4 n) (Z /Xl ...Z Mn ,^,L /x ) = 0. (75) 



The initial condition is again given by the mean-field answer which in d = 2 is 

4 n) (A = l,g R = g )= 5o " f (i Ml . . .£„„) " 2 . (76) 
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The characteristic equations are 

Ts = ' (77) 

d 9R 9r {7r , 

^= (an+^n-ljX), ( 79 ) 

with the boundary conditions as in eq. (|69|) with d replaced by 2. These can again be solved explicitly at s = to 
give: 



i(n-l) /o _ x -i(n-l) 



G n R (A, to) = II / J + ^ log A (AL^.) -2 " (v + lo S A ) " > ( 8 °) 



or, in terms of the rescaled lengths, Li : 



Hn-l) /o„ , \-|(«-D 



The large mass limit corresponds to Li/L^. — > oo in which case 



i=l 



I(n-l) / s \ -i(n-l) 



^...ww-iW^^Ct) ~ r& " ■ (82> 



To recover the asymptotic behaviour in mass space, it is again necessary to take inverse Laplace transforms with 
respect to the fli as fli — > 0. Recalling the definition, eq. of L M we can write 

/n / I - I — j-z — \ n ,~ x |n(n-l) / \ — |n(n— 1) 

JJ e -*«A« dAi ( ^__i ogTO . A .^_^Lj ^__i ogm . A .j - . (83) 

By introducing scaling variables X{ = rhifli and keeping leading order terms in log rhi/rrii the asymptotic behaviour 
of this integral as the rhi — > oo is shown to be 

n / /— = ~ — 3 \ / ~ \ _ i( n_1 ) 

C?(mi...^ > te)~C( ffJi )n[^Iog^mr , J (log^j _ ' (84) 

The density, n = 1, picks up a square root of a logarithm coming from the renormalisation of the reaction rate. However 
the higher order correlation functions pick up additional logarithmic corrections which come from the anomalous 
dimension of the two point function. Note that in d = 2 the asymptotic behaviour retains some memory of the low 
mass cut-offs, ra$. Furthermore the prefactor, denoted above by C{gn) remains dependent on the value of qr, unlike 
in d < 2. 

VIII. RENORMALISED SMOLUCHOWSKI EQUATION 

Although the correspondence between constant kernel aggregation and the A + A — » A system is convenient for 
calculations and greatly simplifies the field theoretic description of the problem, the physics is sometimes obscured 
by this mapping. In this section we show that the results for the average mass density, derived in section IVII CI 
from a field theoretic perspective, can also be interpreted in a more physically transparent way using some heuristic 
arguments closely related to the so-called Smoluchowski approximation used in [23 to take into account fluctuations 
in low dimensional heterogeneous annihilation. 

The idea is as follows. If we are interested in the average mass density, then the only diagrams remaining in 
the field theory associated with the stochastic rate equation, eq. ([Tlj). after all trees have been summed are those 
which renormalise the reaction rate. Therefore it should be possible to understand the problem directly in terms of a 
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standard Smoluchowski kinetic theory with a suitably modified kernel, without reference to the Laplace-transformed 
field R^. Furthermore, the stationary state of this renormalised Smoluchowski theory can be found directly without 
need for any R.G. analysis using a technique, known as the Zakharov Transformation, borrowed from the the theory 
of wave turbulence. 

The variation of the dimensionless reaction rate, g, as the length scale L M is changed is determined by the f3- 
functions computed in section IVII Bl From eq. (ffify and eq. (|KUj) we can obtain the renormalisation law for the 
physical reaction rate, A^(L M ) = g#(L M ) L~ e . Since the Laplace variable, /i, can be thought of as an inverse mass, we 
can make the heuristic substitution, 

L ^{J£) (85) 

to motivate the following mass-dependent reaction rates : 

X R (m) = ^— ^ , d<2, (86) 

XrM = A 7 r , d=2. (87) 



The renormalised Smoluchowski equation (RSE) is then obtained from the mean field Smoluchowski equation by 
substituting the above mass-dependent reaction rate for A. The density should therefore satisfy the following equation 
at large times 



dN(m,t) 

at 



oo 







oo 







poo 

Jo 



dmidm2\R(rn)N(rni,t)N(rri2,t)5(rn — mi — 777-2) 

dmidm2\R(m2)N(m, t)N(mi,t)5(m2 — m — mi) (88) 
dmidm2\R{mi)N(m, i)N(rri2, t)5(m\ — m 2 — rn) 



J r, 

+ — 0(777 — mo). 

777 

The stationary state of this equation is best studied using the method of Zakharov transformations as detailed in 
0,[2(|. For d < 2 and m ^> mo, A^m) = 2ire( J / D)^T^ m~^T^ . The constant flux solution obtained by applying the 
Zakharov transformation to eq. with this kernel is 

d 

/ J \ d+2 2d+2 

N(m)=c K (— J m "3+*" . (89) 

The constant, ck-, can be calculated exactly. However since we have only calculated the renormalised kernel heuris- 
tically, it does not make sense to do so at this point. Note that we recover the predictions of the K41 theory for the 
large mass behaviour of the solution to the RSE. Conversely, for small masses we see that A j r(tt7) = A. In this limit, 
the Zakharov transformation gives the original mean field solution for constant kernel aggregation, 



N(m) = ^^_m-i, (90) 

as found already from our RG analysis. It is an easy task to show that both the mean field and renormalised mean 
field density distributions are local in the sense that the inertial range mass transfer in the stationary state does not 
depend strongly on the source (or the sink which removes large masses if one is present). Therefore both are physically 
realisable solutions and both carry the same flux, J, of mass from small masses to large. The complete stationary 
distribution for constant kernel aggregation should therefore exhibit a crossover from the mean field solution to the 

renormalised mean field solution at a mass, tt7 c which is given by {J / D)~^m c d+2 «1. We note that this crossover 
from a mean field spectrum to a fluctuation dominated one has been conjectured to occur in other turbulent systems, in 
particular in wave turbulence [27j| . although these systems do not readily lend themselves to such systematic analysis. 
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In d = 2, analysis of the RSE equation allows us to obtain the logarithmic correction to the mean field s pect rum, 
again without resort to R.G. arguments. Strictly speaking, the Zakharov transformation technology of [Tol l26j only 
works for homogeneous kernels which leaves us with the question of what to do with the logarithm in eq. (|87[). It 
turns out that the approach is easily adapted to extract the logarithmic correction to the spectrum at large masses. 
We give a brief outline here. In d — 2, the renormalised interaction, eq. (|87[). behaves for large masses as 



X R (m) 



1 



In 



\m J 



(91) 



After application of the Zakharov transformations described in 10], the stationary RSE, eq. (|88|) . it becomes 



f 

Jo 



dm\drri2 



A—) 

\ni2rn0 j 



In 



\ _1 



^rairao 
S(m — mi — 7712). 



ln( — ) Mmi)7V(m 2 ) 
\m J 

N(m) iV(^— !-) fllV 
m 2 \rri2 J 

1 / x 2^ 

N(m)N( -) — 

mi \mi 



(92) 



Let us now look for a solution of the form N(m) = ck \n(m/mo) y m~ x . This substitution yields a rather messy 
expression which we analyse by introducing new integration variables, /ii, \±2 defined by vrt\ = m/ii, m 2 = m/i 2 and 
then expanding the resulting expression as a power series in ln(m/mo) _1 . After some algebra one obtains : 
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(93) 



It is clear that the leading logarithms cancel out asymptotically as m — » 00 if we choose y — 1/2. The integrand then 
vanishes asymptotically for x = 3/2 as in the usual mean field case. The renormalised Kolmogorov spectrum in d = 2 
for large masses is therefore 



P(m) c K 




(94) 



as found from the RG analysis of section IVlI Fl 

Let us close this section by discussing the connection between the renormalised Smoluchowski equation and the 
Smoluchowski approximation used in [25( to study the kinetics of heterogeneous annihilation. The essence of the 
argument used in 25] is as follows. Consider a heterogeneous system of annihilating particles where the reactants 
have a continuous distribution of diffusivities with the slower particles being less probable than the faster ones. Now 
consider the reaction between particles with diffusivity D and "slower" particles with diffusivity D 1 < D. Since 
the slower particles are rare we can neglect reactions between slow particles. We therefore estimate the effective 
reaction rate for the slow particles by considering each slow particle to be stationary in a uniform background cloud 
of faster particles and calculating the diffusive flux of fast particles reaching the slow particle from the background 
cloud. One estimates the effective reaction rate in d dimensions at time t to be (D + ]j'Y^ 2 t~ 1Jrd ^ 2 . Substitution 
of this effective reaction rate in the mean field rate equation constitutes the Smoluchowski approximation which 
greatly improves the estimate of the asymptotic decay rate of the particle density in d < 2. In the case at hand, 
the dynamics naturally generates a distribution of particle masses in which heavy particles are much rarer than the 
lighter ones. So, although in our model all particles have equal diffusivity, D, in the reference frame of a heavy 
particle we can consider the main interaction to be with a uniform background cloud of light particles. Thus one 
can envisage an effective reaction rate for the large mass particles of (2L>) d / 2 £- 1+d / 2 . However in our case the mass 
flux into the system is a constant, J, so that the system reaches a stationary state in which the natural unit of time 
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for particles of mass m is t ~ (J~ 1 D~ d ^ 2 m) 2 ^ d+2 \ Thus the time dependent effective reaction rate of [2^ should 
be replaced in the stationary state of the constant kernel aggregation problem with a mass dependent reaction rate 
A j R(m) ~ ]J d / 2 (J]J d / 2 }-( d - 2 )/( d + 2 ) m ( d - 2 )/( d + 2 ) which we recognise as the large mass behaviour of the renormalised 
kernel, eq. Thus the large mass asymptotics of the RSE corresponds to the correction of the mean field rate 

equation by the Smoluchowski approximation. The field theory approach with which we derived the original results 
allows unambiguous identification of the set of diagrams which have been summed to give this approximation. 



IX. CONCLUSIONS 



To summarise, the most important result of this work is the following fact : the steady state mass density PDF of 
cluster-cluster aggregation in the presence of a steady source of monomers exhibits non-trivial multiscaling in dimen- 
sions less than or equal to 2. Technically, this means that the exponents, £ n , describing the large mass asymptotics 
of the density correlation functions, C n (m\ . . .m n ) = (N mi . ..iV mn ), depend nonlinearly on the order, n, of the 
correlation function. Physically it means that the system is in a regime where strong correlations (or in this case, 
anticorrelations) between particles dominate the statistics putting the problem outside of the domain of applicability 
of mean-field theory. In the early sections of the article we established the presence of multiscaling more-or-less 
rigorously by calculating the asymptotic scaling behaviour of the C n for n = 0, 1, 2 and remarking that they do not 
lie on one line. The latter part of the article is devoted to detailing a nonrigorous derivation of ( n for general n as an 
expansion in e = 2 — d using renormalisation group techniques. To leading order in e, the scaling curve is quadratic 
in n. In d = 2 the power law corrections to the mean-field exponents become logarithmic. Our results are interesting 
both from the perspective of particle systems and from the perspective of turbulence. In the context of aggregation, 
the understanding of the detailed role of fluctuations beyond their effect on the density is still embryonic. The effect of 
fluctuations on the density can be taken into account using the so-called Smoluchowski approximation. In the closing 
section of the article we showed how our approach is connected to this approximation. Whereas the Smoluchowski 
approximation is only useful for the density, our methods work for arbitrary correlation functions. The connection to 
the theory of turbulence stems from the analogy between the mass cascade in the stationary state of cluster-cluster 
aggregation and the energy cascade in the stationary state of homogeneous isotropic turbulence. Throughout the 
article we have tried to emphasise both the usefulness and limitations of this analogy. Both systems exhibit nontrivial 
multiscaling. In the case of turbulence, this multiscaling has proven to be very difficult to understand quantitatively. 
One might hope that the analogy developed here might, if viewed in the correct way, provide some insight. At present, 
however, this is merely an aspiration. 
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APPENDIX A: EFFECTIVE HAMILTONIAN OF MASS MODEL. 

A microstate of the mass model is given by the vector of occupation numbers, N = {^V^ m }^ GZ d ?mGZ +. The 
probability measure Pt(iV) on the space of microstates satisfies the Master Equation, which follows from the dynamics 
of the mass model: 

§i p t(ti) = Y d E E ((Ns,m + l)Pt({. ..ns, m + 1. . .n*/, m - 1. . .}) - N^ m P t {N)) 

m (x,x') 

+ ^E E ({N^ mi + l){N^ m2 + l)P t {{...N^ mi +l...N Stm2 +l...N^ m -l...}) - N Stmi N^ m2 P t (N)^ 

x ra=rai+ra2 

E E (—8(m - m )Pt ({• • • N s<m - 1 . . .}) - — S(m - m )Pt(N)) ■ 

x rn x ' 

The first line of the Master Equation describes diffusion, second - aggregation, the third - deposition. By " . . ." 
we denote components of the microstate, which are the same on the right hand side as on the left hand side. Let 
\P t ) — Pt(N) Ylg m (at rn ) Ns,rn |0) be the state vector corresponding to probability measure P t . 
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The evolution equation for \P t ) follows by differentiating \P t ) with respect to time and using the master equation 
to express the result back in terms of \P t ). A simple calculation shows that 

d t \P t ) = -H tm \P t ), (Al) 

where H m is an effective Hamiltonian of the mass model given by : 

H tm = D Yl (4,m - a S',m)( a S,m ~ a S > , ro ) 

-A ^ ^2 ^ m ~ mi ~ m 2) a £ i m a $,™>i a $,™>2 



x m,m 1 ,rri2 



-E^^-^)(4,^-l) (A2) 

where (5, x') denotes summation over pairs of nearest neighbours. Note that we replaced A with 2A in order to 
simplify the notation. The Hamiltonian operator, eq. (|A2[) . already includes the "shift" introduced to simplify the 
computation of correlation functions as described in appendix 151 



APPENDIX B: RELATION BETWEEN CORRELATION FUNCTIONS OF INTERACTING PARTICLE 

SYSTEMS AND STOCHASTIC RATE EQUATIONS. 

In this appendix, we establish the relation between interacting particle systems and stochastic rate equations. To 
avoid unnecessary complications, we start with an example of a zero-dimensional system. We then show how the zero 
dimensional results get modified in d > using the example of the mass model. 

The microstate of a zero dimensional particle system is determined by the number of particles present at time 
t. Let the probability that there are TV particles left at time t be P t (N). If the particle dynamics is Markovian, 
this probability satisfies the linear master equation which is first order in time. The principle step of Doi's map 
is to rewrite the master equation as an evolution equation for a certain state vector in the Fock space spanned by 
microstates of the particle system. The Fock space T can be introduced as follows. Let \N) be an element (a "state") 
of JF, corresponding to a microstate with TV particles, Af = 0,l,2, — The state |0) is often referred to as 'vacuum'. 
Then T is defined as a space of all linear combinations of states \N) with complex coefficients. Let a) be the creation 
operator acting in T . By definition, 

a^\N) = \N + 1) (Bl) 

The annihilation operator, a, acting in T is defined as follows: 

[a, at] = 1, (B2) 
a\0) = 0, (B3) 

where [A, B] = AB — BA denotes the commutator of the operators A, B. It then follows that 

a\N) =N\N-1). (B4) 

It is also convenient to define the left vacuum (0| (an element of the space dual to J 7 ) by means of the following 
relations: 

(0|0) = 1, 

(0|a f = 0. (B5) 

Consider the following element of T\ 



oo 

\P t ) = J2 Pt(N)\N). 

N=0 



(B6) 
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The state \P t ) can be used to generate all probabilities P t (N), as Pt(N) = (0\a N \P t ). The vector \P t ) satisfies an 
equation, which follows from the master equation for probabilities: 

j t \P t ) = -H\P t ), (B7) 

where H is an evolution operator (Hamiltonian). The exact form of the Hamiltonian follows from the microscopic 
rules of evolution. For example, if particles undergo pairwise annihilations at rate A, then (see, for example, jl8j |) 

= -A(a 2 -a+ V). (B8) 

The first term on the right hand side of (|B8|) describes annihilation, while the second term comes from the "minus" 
term in the master equation, which accounts for the probability of non-reaction. The formal solution to (|B7|) is 

\P t ) = e~ tH \P ), (B9) 

where the state \Pq) is determined by the initial distribution of particles. The problem of solving Master Equation is 
now reduced to solving an effective Schroedinger Equation in imaginary time, which can be often done using powerful 
methods of quantum/statistical field theory. 

Another important observation of Doi's theory is the averaging formula. Let Z t (J) be the generating function 
of moments of the probability distribution P t (N). By definition, Z t (J) = E(e JNt ) = e JN P t (N). Using the 
definitions above, one can prove the following: 

Z t (J) = (0\e a e Jafa \P t ). (BIO) 
To show this, expand the right hand side to give 

Zt(J) = E E Pt{N) J —{0\e a ^a) K \N) 

N=0 K=0 ' 



OO OO yjfc 

E E j^N K P t (N)(0\e a \N) 



N=0K=0 

oo oo 



= E E J ^ K ?m 



K\ 

7V=0 K=0 



oo 



Y,e JN P t (N) = E(e JN ). 



N=0 



In the proof we used the fact that TV = a) a is the operator measuring the number of particles, N\N) = N\N), and 
the identity e a a^ N = (a) + l) N e a . The form of the correlation function on the right hand side of eq. (|B1Q|) is not very 
suitable for practical computations due to the presence of the shift operator e a inside the brackets. This problem can 
be overcome by commuting e a to the right and using e a |0) = |0). It follows from the commutation relations [eq. (|B2|) ] 
that e a O(a\ a) = 0{a^ + 1, a) for any operator 0(a) , a). Using this fact one finds that 

Z t (J) = {0\e J ^^ a e- t& \P ), (Bll) 

where H(a^a) = H(a) + l,a), \Pq) = e a |Po)- The expression ^o|e J ( at+a ) a can be simplified further. Note that 
[a) a, a] = —a. In other words, the operators a and a) a form a basis of a Lie algebra isomorphic to a subalgebra of 
s/(2) consisting of upper triangular matrices. Consequently, the Campbell-Hausdorff formula implies that 

e J(a^a)a = e Ja^a e f(J)a^ ( B12 ) 

where f(J) is a function to be determined. Differentiating both sides of eq. (|B12|) with respect to J, commuting all 
operators multiplying exponents in the derivatives to the right, and comparing both sides of the resulting equality, 
we find a differential equation for /: 



f(j) = f(J) + 1, 



(B13) 
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which should be solved with the boundary condition /(0) = 0. The answer is f(J) = e J — 1. Taking into account 
that (0|e Jata = (0|, we find 

Z t (J) = (Ole^-^e-^lPo), (B14) 

The right hand side of this relation can be rewritten as a path integral using the Trotter formula. Assume for simplicity 
that initial probability distribution of the number of particles Po(t) is Poisson with intensity No. Then, 

where the integration is performed over the space of complex paths (a(r), a(r)) re [ ,t] an d 

S eff {t) = J^dr ^d(r)d T a(r) + H(a(r), a(r)) - A^(r)) (B16) 

is an effective 'action functional'. Here H (a(r) , a(r)) is the symbol of operator H(a\a), which is assumed to be 
normally ordered. We conclude that 

E(e JN >) Pt = < e (''-iW*)> Se//(t)) (B17) 

where we stressed by our notations, that the averaging in the left hand side of eq. (|B17|) is performed over the space 
of microstates using probability distribution P t , while averaging in the right hand side is performed over the space of 
all paths (a(r), a(r)) rG [o j t] using the functional measure e~ Sef f^\ Differentiating both sides of eq. (|B17j) with respect 
to 6 = e J — 1 and setting = we find that moments of a(t) correspond to factorial moments of P t : 

E(N t (N t - l)(N t - 2)...(N t —k + l)) Pt = (a(t) k ) Seff{t) ,k = 1, 2, 3, . . . (B18) 

In part icular E(N t ) = (a(t)),E(N 2 ) = (a(t) 2 + a(t)), etc. 

Eq. (|B18|) allows one to capture strong non-mean field behaviour of an interacting particle system even if the field 
theory characterised by eq. (|B16|) is well approximated by mean field theory. The simplest example is as follows. 
Assume that (a(t) k ) « (a(t)) k , but E(N t ) = (a(t)) <C 1. Using these assumptions and eq. (|B18|) to evaluate moments 
of Nt, we find that E(N k ) = E(N), which is essentially non-mean field behaviour. 

For particle systems with pairwise local interactions, eqs. (|B17|) and (|B18|) can be formulated in terms of stochastic 
differential equations thus avoiding references to non-rigorous path integral constructions. To illustrate this, consider 
the reaction A + A — » in zero dimensions. The effective Hamiltonian is given by eq. dBEJ), H = H(a^ + l,a) = 
\{2a^a 2 + a) 2 a 2 ). The corresponding functional integral measure can be rewritten using the Hubbard- Stratonovich 
identity as 

- P dr( a(r)^a(r)+A(2ata 2 +at 2 a 2 )-7V a(r) J , ft 2 

- Pdr( a(r)a T a(r)+2Aa t o 2 -Aroo(r)-i72A^(r)a t (r)a(r) ] 

x e ° V /. (B19) 

The field ^(t) is standard Gaussian white noise. Note that the expression in the exponent in the second line of 
eq. (|B19|) is linear in the field a\ which allows it to be integrated out. This results in a (5-functional with an argument 

d T a(r) + 2Xa 2 - N S(r) - iV\^(r)a(r) = 0, (B20) 

in which we recognise a rate equation of A + A — > augmented by imaginary multiplicative noise term (Lee-Cardy 
equation). 

We can now interpret eqs. (|B17|) and (|B18|) as follows. Let Nt be the number of particles left after time t in the 
system of annihilating particles. Let a(t) be the solution to stochastic differential eq. (|B20|) . Then eqs. (|B17|) and (|B18|) 
are true, given that (...)s eff denotes averaging over the white noise £(£) |2^. In the stated form, the correspondence 
between the Markov chain describing A + A — > system and the Lee-Cardy stochastic partial differential equation 
can be proven rigorously and extended to a large class of interacting particle systems in d > [2^ . 
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Here we will only give a path integral derivation of the relation between correlation functions in the mass model in 
d > and correlation functions of Stochastic Smoluchowski equation, eq. 

First, consider the case of discrete space and mass. A microstate of the mass model is specified by stating the 
number of particles of a given mass at a given site. In other words, the microstate is a vector {N$ irn }£ e z d ,mez + with 
non- negative integer components. The operator of mass distribution is 

^x,m ^jt rn ^ j x,rn 5 (B21) 

where at and a^ m are creation and annihilation operators satisfying the following commutation relations: 

[4,m' a x',m'] = 5x,x'5m,m'- (B22) 

Similarly to the zero dimensional case considered above, the generating functional for mass distribution correlation 
functions can be written as follows: 

Z t [j\ EE E^j^^rr,^^^^ 

= (ole^-^e^-^^^e-^^lPo), (B23) 

where 

Hum = (4,m ~ 4',m)( a *,™ - a*',m) 

m (x,x') 

XV 171,1711,1712 



+ 4,mi4,m 2 a ^™i a ^' 



X mi ,1712 



-E^- <5 ( m - m °)(4,, m -l) (B24) 

x,m 

is an effective Hamiltonian of the mass model and Po is an initial probability measure on the space of microstates. 
An outline of the derivation of Hmm is presented in appendix 1X1 

As before, we need to commute the projection operator exp^^ rn a x,m to the right. Computations, which are 
completely analogous to those performed in O-dimensional case give: 

Z t [J] = (0|e^- (exp(J ^ m) " 1)a ^' m e- #MM |Po), (B25) 

where Hmm is obtained from Hmm by the shift at m — > at m + 1 and |Po) = > m ax,m |Po). The right hand side of 
eq. (|B25|) can be re-written in path integral form in the conventional way. A lengthy but straightforward calculation 
shows that 



t [j) = / n 

x'm'r' 



d(j)(x\ m', r r )d(j)(x r , m', r) 



e T,z,J ex P( J ^™)-±)Hx,rn,t) e _s M M(t)^ ^ B26 ^ 



where 



5mm (*) = J dt [Yl 

^ \x,17l 



d 

— -DA)</>(x,m,t) 



<j>{x, m, t)d t ct)(x, m, t) + H M m[(/>, <t>] 

J dt^ \ 4>{x, m, t) 
—A ra', i)(/)(x, m — m', t) 

+2A0(f , m, t) V" </>(f , ra', t) - — 5 m , mo 

mo 

m' 

+A 0(5, ra, t)4>(x, ra, t)(j>(x, ra/, t)4>(x, ra/, t) j , (B27) 
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where A is a discrete Laplacian and J is the rate of input of mass into the system. Note that the expression in square 
brackets in the right hand side of eq. (|B27j) is just the constant kernel Smoluchowski equation. The last term in 5mm 
accounts for all correlation effects. The exponential of this term can be rewritten using the Hubbard- Stratonovich 
transformation as follows: 



j \{di{x'y)e-^l dt ^^ t)2 

x' ,t' 



(B28) 



Note that the field £ is Gaussian, uncorrelated both in space and time. Using eq. (|B28[) . the functional measure of 
integration in eq. (|B26|) can be rewritten in the form 



J II dttfy^-lti d ^^ e -Io^, m ^ m ^ L ^\ (B29) 



where 



L[</>,£] = (-Q t -DA\(t>(m)-\ <t>{x,m',t)<i>{x,m-m',t) 



+2A0(f , m, t) V (j)(x , m', t) - —S m mo - iV2\(j)(m)r](x, t). (B30) 

The exponent in the right hand side of eq. (|B29|) is linear in 0. Hence the path integral over fields 0, and £ localises 
to paths satisfying Euler-Lagrange equation 



/ dtV^>,t)L[«], (B31) 
7o -* _ 



or 

L[0,£]=O, (B32) 

which is a discrete version of the Smoluchowski equation, eq. JTT}. For the sake of clarity, we restate the result 
concerning the relation between the mass model and stochastic Smoluchowski equation here. In order to calculate the 
generating functional of density correlation functions in the mass model [eq. (|B23|) ]. one has to solve the Stochastic 

Smoluchowski equation, eo. (|B32|) . for m, t), then average e^ s > m ^ exp ^ Jx,m ^ 1 ^^ x,m ^ with respect to Gaussian 

white noise £. In other words, 

Z t (J) = s^ e ^.- (CTp(Jfif ' m) - 1)0K](5r,m) ^ (B33) 

Our final task is to discuss the modification and consequences of eq. (|B33|) in the continuous limit. The latter is 
taken according to the following set of rules: 

x 

x -> — , 
a x 
m 

m — > — , 

Clm 

. , x m 1 _ v 
0( — , — ,r) -> 0(x,m,t), 

d x d m 
n n n a n 



Qj X CLm, d^d* 

D 



x {Jj rn Lt x ujm 

D 



J4 
J -> — 

(2 777, 

m - ^, (B34) 

CI777 
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where a x and a m are lattice cut-offs in x- and ra-spaces respectively. The continuous limit is obtained by performing 
replacements eq. (|B34|) in eqs. (|B30|) and (|B33|) and taking lattice cut-offs to zero while keeping other parameters 
fixed. As a result one recovers the stochastic Smoluchowski equation (SSE) [eq. JUJ]. 

Note that the continuous field theory equivalent to SSE is renormalizable in dimensions two and less, therefore 
eq. (|B34|) is justified in these dimensions only. The continuous counterpart of eq. (|B33|) is 

Z t [J) =E^e<f / d ^ m ( e ^( J (^ m ))- 1 )^ ) [^(^' m )^ (B35) 



E\4(x,m,t)J = E^N t (x,m)J, 
E[4>(xi,mi,t)(t>(x2,m2,t)) = E[ N t (x u mi)N t (x 2 , m 2 ) 



APPENDIX C: EXPRESSION OF PROBABILITY OF MULTI-PARTICLE CONFIGURATIONS IN 

TERMS OF SOLUTIONS TO SSE. 

In this appendix, we show how probability of multiparticle configurations can be calculated from SSE. Equa- 
tion (|B35|) leads to the following relation between correlation functions: 



-S d (x 1 -x 2 )S{rm -m 2 )Nt(x 1 ,rn 1 ) S j , (CI) 

and so on. Suppose that we are interested in the statistics of the total number of particles AN t (x,m) in a volume 
element AV centred around x with masses in the interval [m, m + Am]. In terms of the local mass distribution, 

p rm-\-5m 

AN t (x,m) = d d x' dm'N t (x',m'). (C2) 

J AV Jrn 

Let 

r rm-\-5m 

A0 t (f,m) = / d d x' dm'(j) t (x' ,m'). (C3) 

J AV Jrn 

Integrating eq. (|C1|) with respect to mass and space, we find 

J5^A0(£,ra,t)^ = E^AN t (x,m,t) 

£^A</> 2 (£,m,f)^ = E^AN t (x,m,t)(AN(x,t,m)-l) 

E(A(j) n (x,m,t)\ = e( Y[(AN(x : t,m)-k)\ (C4) 
^ ' ^k=o ' 

which is a multi-dimensional counterpart of O-dimensional result of eq. (|B18[) . 

In this paper we study scaling properties of probability of finding multiple particles of large mass in A V Am. Density 
of such particles is low. Thus, factorial moments entering the right hand side of eq. (|C4|) can be estimated in the limit 
of large mass m as follows: 

✓ ri— 1 \ oo / n— 1 \ / \ 

e( Y[(AN(x,t,m) - k) J ^ ( JJ(p — fc) )Prob(xi, ...,x p G AV m ,mi,...,m p G [m,m + Am] 

^ k=0 ' p=n ^ k=0 ' ^ ' 

w nlProb fi,...,f n G AV ; mi , . . . , m p G [m, m + Am] ) (C5) 



Combining this result with eq. (|C4[) . we obtain desired relation between probabilities of multi-particle configurations 
and moments of solutions to SSE in the limit of large masses: 



Probl xi, . . . , x n G AV; mi, . . . , m v G [m, m + Am] ) = —E [ A(j) n (x, m, t) 



(C6) 
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